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Abstract. The celebrated Kuramoto model captures various synchronization phenomena in 
biological and man-made dynamical systems of coupled oscillators. It is well-known that there exists 
a critical coupling strength among the oscillators at which a phase transition from incoherency to 
synchronization occurs. This paper features four contributions. First, we characterize and distinguish 
the different notions of synchronization used throughout the literature and formally introduce the 
concept of phase cohesiveness as an analysis tool and performance index for synchronization. Second, 
we review the vast literature providing necessary, sufficient, implicit, and explicit estimates of the 
critical coupling strength in the finite and infinite-dimensional case and for both first-order and 
second-order Kuramoto models. Third, we present the first explicit necessary and sufficient condition 
on the critical coupling strength to achieve synchronization in the finite-dimensional Kuramoto model 
for an arbitrary distribution of the natural frequencies. The multiplicative gap in the synchronization 
condition yields a practical stability result determining the admissible initial and the guaranteed 
ultimate phase cohesiveness as well as the guaranteed asymptotic magnitude of the order parameter. 
As supplementary results, we provide a statistical comparison of our synchronization condition with 
other conditions proposed in the literature, and we show that our results also hold for switching and 
smoothly time-varying natural frequencies. Fourth and finally, we extend our analysis to multi-rate 
Kuramoto models consisting of second-order Kuramoto oscillators with inertia and viscous damping 
together with first-order Kuramoto oscillators with multiple time constants. We prove that such a 
heterogenous network is locally topologically conjugate to a first-order Kuramoto model with scaled 
natural frequencies. Finally, we present necessary and sufficient conditions for almost global phase 
synchronization and local frequency synchronization in the multi-rate Kuramoto model. Interestingly, 
our provably correct synchronization conditions do not depend on the inertiae which contradicts prior 
observations on the role of inertial effects in synchronization of second-order Kuramoto oscillators. 
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1. Introduction. A classic and celebrated model for the synchronization of cou- 
pled oscillators is due to Yoshiki Kuramoto [33] . The Kuramoto model considers n> 2 
coupled oscillators each represented by a phase variable 9i G T^, the 1-tours, and a 
natural frequency uji G M. The system of coupled oscillators obeys the dynamics 

K " 

e^^uj,- -Y.^m{9,~9,), ie{l,...,n}, (1.1) 

" i=i 

where K > is the coupling strength among the oscillators. 



The Kuramoto model (1.1) finds application in various biological synchroniza- 
tion phenomena, and we refer the reader to the excellent reviews [HI [2] for various 
references. Recent technological applications of the Kuramoto model include motion 
coordination of particles [13] i synchronization in coupled Josephson junctions |5lj, 
transient stability analysis of power networks [20], and deep brain stimulation [47] . 

The Critical Coupling Strength. Yoshiki Kuramoto himself analyzed the 



model (1.1) based on the order parameter re^^ = ^ X]j=i ^'^^j which corresponds the 
centroid of all oscillators when represented as points on the unit circle in C^. The 
magnitude of the order parameter can be understood as a measure of synchronization. 
If the angles Oi{t) of all oscillators are identical, then r = 1, and if all oscillators are 
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spaced equally on the unit circle (splay state), then r = 0. With the help of the order 



parameter, the Kuramoto model (1.1 1 can be written in the insightful form 



^oj,^ Krsm{0,~ij), te{l,...,n}. (1.2) 



Equation (1.2) gives the intuition that the oscillators synchronize by coupling to a 
mean field represented by the order parameter re"''. Intuitively, for small coupling 
strength K each oscillator rotates with its natural frequency Wi, whereas for large 
coupling strength K all angles 9i{t) will be entrained by the mean field re"'' and the 
oscillators synchronize. The threshold from incoherency to synchronization occurs 
for some critical coupling i^Tcriticai- This phase transition has been the source of 
numerous research papers starting with Kuramoto's own insightful and ingenuous 



analysis [331l31j. For instance, since r < 1, no solution of (1.2 1 of the form di(t) = 
dj{t) can exist ii K < \uJi — ujj\/2. Hence, K > \uJi — Wj|/2 provides a necessary 
synchronization condition and a lower bound for ii^criticai- Various various necessary, 
sufficient, implicit, and explicit estimates of the critical coupling strength -/^critical for 
both the on-set as well as the ultimate stage of synchronization have been derived in 
the vast literature on the Kuramoto model [SlillliaiSnilSllSllSTlllHlIIllMlIini 
SlIIllinilMlllllSllilSnilHllMlllSl-To date , only explicit and sufficient (or necessary) 
bounds are known for the critical coupling strength i^criticai h^ the Kuramoto model 
(1.1), and implicit formulae are available to compute the exact value of i^Tcriticai- 



The Multi-Rate Kuramoto Model. As second relevant coupled-oscillator 
model, consider m > second-order Kuramoto oscillators with inertia and viscous 
damping and n — m > first-order Kuramoto oscillators with multiple time con- 
stants. The multi-rate Kuramoto model evolving on T" x M™ then reads as 



K " 

MiO, + DA = sm{e, - e/j) , ie {l,...,m}, 



n 
1=1 



K " 

DiOi = uji sm{9i — 9j) , i E {m + 1, . . . , n} . 



(1.3) 



n 
1=1 



where Mi > 0, Di > 0, and cj^ e M for i e {1, . . . , n} and K > 0. Note that we allow 



for m G {0, n} such that the model ( 1.3 ) is of purely first or second order, respectively. 

The multi-rate Kuramoto model ( |1.3[ ) finds explicit application in the classic 
structure-preserving power network model proposed in [7]. For m — n, the model 



(1.3) is a purely second-order system of coupled, damped, and driven pendula, which 
has been used, for example, to model synchronization in a population of fireflies |23j . 
in coupled Josephson junctions [5T], and in network-reduced power system models [TT] . 
For m = n, unit damping Di = 1, and uniform inertia Mi = M > 0, the 



second-order Kuramoto system (1.3) has been extensively studied in the literature 
[T51 H51 [551 123 HI H] • The cited results on the inertial effects on synchronization 
are controversial and report that synchronization is either enhanced or inhibited by 
sufficiently large (or also sufficiently small) inertia M. For the general multi-rate 



Kuramoto model (1.3) no exact synchronization conditions are known. 



1.1. Contributions. The contributions of this paper are four-fold. First, we 
characterize, distinguish, and relate different concepts of synchronization and their 
analysis methods, which are studied and employed in the networked control, physics, 
and dynamical systems communities. In particular, we review the concepts of phase 
synchronization and frequency synchronization, and introduce the notion of phase 
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cohesiveness. In essence, a solution to the Kuramoto model (1.1) is phase cohesive if 
all angles are bounded within a (possibly rotating) arc of fixed length. The notion of 
phase cohesiveness provides a powerful analysis tool for synchronization and can be 
understood as a performance index for synchronization similar to the order parameter. 

As second contribution, we review the extensive literature on the Kuramoto 
model, and present various necessary, sufficient, implicit, and explicit estimates of 
the critical coupling strength for the finite and infinite-dimensional Kuramoto model 
in a unified language [S3lll[ail3[2nilMl[2aEZlllSl[Il[2a[iaill[Illia[S 
El [35l m [50l [8j [24l [26] . Aside from the comparison of the different estimates of the 
critical coupling strength, the second purpose of this review is the comparison of the 
different analysis techniques. Furthermore, we briefly survey the controversial results 
[131 [ini [13 [551 [23 [3 [2 on the role of inertia in second-order Kuramoto models. 

As third contribution of this paper, we provide an explicit necessary and suffi- 
cient condition on the critical coupling strength to achieve exponential synchroniza- 
tion in the finite-dimensional Kuramoto model for an arbitrary distribution of the 



natural frequencies w^, see Theorem 4.1 In particular, synchronization occurs for 
K > iiTcriticai = "^max " i^min, whcrc Wmax and Wmin are the maximum and minimum 
natural frequency, respectively. The multiplicative gap ii'criticai/-?^ determines the 
admissible initial and the guaranteed ultimate level of phase cohesiveness as well as 
the guaranteed asymptotic magnitude r of the order parameter. In particular, the 
ultimate level of phase cohesiveness can be made arbitrary small by increasing the mul- 
tiplicative gap Kcnticax/ K . This result resembles the concept of practical stability in 
dynamics and control if K and -K^criticai are understood as a synchronization-enhancing 
gain and as a measure for the desynchronizing non-uniformity among the oscillators. 
Additionally, our main result includes estimates on the exponential synchronization 
rate for phase and frequency synchronization. We further provide two supplementary 
results on our synchronization condition. In statistical studies, we compare our con- 
dition to other necessary and explicit or implicit and exact conditions proposed in 
the literature. Finally, we show that our analysis and the resulting synchronization 
conditions also hold for switching and smoothly time-varying natural frequencies. 

As fourth and final contribution, we extend our main result on the classic Ku- 
ramoto model (1.1) to the multi-rate Kuramoto model (1.3). We prove a general 
result that relates the equilibria and local stability properties of forced gradient-like 
systems to those of dissipative Hamiltonian systems together with gradient-like dy- 



namics and external forcing, see Theorem 5.1 As special case, we are able to show 



that the multi-rate Kuramoto model is locally topologically conjugate to a first-order 
Kuramoto model with scaled natural frequencies, see Theorem |5.3| Finally, we present 
necessary and sufficient conditions for almost global stability of phase synchronization 
and local stability of frequency synchronization in the multi-rate Kuramoto model, 
see Theorem |5.5[ Interestingly, the inertial coefficients Mi do not affect the syn- 
chronization conditions and the asymptotic synchronization frequency. Moreover, the 
location and local stability properties of all equilibria are independent of the inertial 
coefficients Mi , and so are all local bifurcations and the the asymptotic magnitude of 
the order parameter. Rather, these quantities depend on the viscous damping param- 
eters Di and the natural frequencies uji. Of course, the inertial terms still affect the 
transient synchronization behavior which lies outside the scope of our analysis. These 
interesting and provably correct findings contradict prior observations on the role of 
inertia inhibiting or enhancing synchronization in second-order Kuramoto models. 

The remainder of this paper is organized as follows. Section [2] reviews different 
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concepts of synchronization and provides a motivating example. Section [3] reviews 
the hterature on the critical coupling strength in the Kuramoto model. Section [4] 
presents a novel, tight, and explicit bound on the critical coupling as well as various 
related properties, performance estimates, statistical studies, and extensions to time- 
varying natural frequencies. Section [5] extends some of these results to the multi-rate 
Kuramoto model. Finally, Section |6jconcludes the paper. 

Notation. The torus is the set = ] — tt, +tt], where —it and +tt are associated 
with each other, an angle is a point G T^, and an arc is a connected subset of T^. 
The product set T" is the n-dimensional torus. With slight abuse of notation, let 
1^1 — ^2! denote the geodesic distance between two angles 9i € and 62 ^ T^- For 
7 € [0, tt], let A(7) C T" be the set of angle arrays . . . , 6n) with the property that 
there exists an arc of length 7 containing all ^i, ... , 6'„ in its interior. Thus, an angle 
array 9 G A(7) satisfies max^ jg{i_...^„} \9i — 9j\ < 7. For 7 £ [0,7r], we also define A(7) 
to be the union of the phase-synchronized set {6' e T" | 9i — 9j, i,j G {1, . . . , n}} and 
the closure of the open set A(7). Hence, 9 G A(7) satisfies max^^^g^i^ „} \9i — 9j \ < 7; 
the case 9 G A(0) corresponds simply to 9 taking value in the phase-synchronized set. 

Given an n-tuple {xi, . . . , Xn), let a; G K" be the associated vector, let diag(a;i) G 
M" be the associated diagonal matrix, and let Xmax and a;min be the maximum and 
minimum elements. The inertia of a matrix A G M"^" are given by the triple 
{i/g, fc, i^u}, where Vs (respectively i^^) denotes the number of stable (respectively 
unstable) eigenvalues of A in the open left (respectively right) complex half plane, 
and Vc denotes the number of center eigenvalues with zero real part. The notation 
blkdiag( Ai , . . . , An ) denotes the block-diagonal matrix with matrix blocks , . . . , An . 
Finally, let /„ be the n-dimensional identity matrix, and let Ipxq and Opxq denote 
the p X q dimensional matrix with unit entries and zero entries, respectively. 

2. Phase Synchronization, Phase Cohesiveness, and Frequency En- 
trainment. Different levels of synchronization are typically distinguished for the 



Kuramoto model (1.1 1. The case when all angles 9i(t) converge exponentially to a 
common angle ^oo G as i — )• 00 is referred to as exponential phase synchronization 
and can only occur if all natural frequencies are identical. If the natural frequencies 
are non-identical, then each pairwise distance \9i{t) — 9j{t)\ can converge to a con- 
stant value, but this value is not necessarily zero. The following concept of phase 
cohesiveness addresses exactly this point. A solution 9 : M>o T" to the Kuramoto 



model (1.1) is phase cohesive if there exists a length 7 G [0,7r[ such that 9{t) G A(7) 
for all t > Q, i.e., at each time t there exists an arc of length 7 containing all angles 
9i{t). A solution 9 : M>o — > T" achieves exponential frequency synchronization if 
all frequencies 9i(t) converge exponentially fast to a common frequency 0oo G K as 
t — >■ CXI. Finally, a solution 9 : M>o — >■ T" achieves exponential synchronization if it is 
phase cohesive and it achieves exponential frequency synchronization. 

If a solution 9{t) achieves exponential frequency synchronization, all phases asymp- 
totically become constant in a rotating coordinate frame with frequency 9oo , or equiv- 
alently, all phase distances \9i{t) — 9j{t)\ asymptotically become constant. Hence, the 
terminology phase locking is sometimes also used in the literature to define a solution 
9 : M>o T" that satisfies 9i(t) = 9^0 for alH G {1, . . . , n} and for alH > |35','35','53] 
or 9,it) - 9j{t) = constant for all i,j G {1, . . . ,ri} and for alH > HmH, 49, 50 . 
Other commonly used terms in the vast synchronization literature include full, ex- 
act, or perfect synchronization for phase synchronizatiorj^ and frequency locking, fre- 



Notc that 2 understands phase locking synonymous to phase synchronization as defined above. 
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quency entrainnient, or partial synchronization for frequency synchronization. 

In the networked control community, boundedness of angular distances and con- 
sensus arguments are typically combined to establish frequency synchronization [141 
[29, J2j 20, 24, 26 . Our latter analysis in Section [4] makes this approach explicit by 
distinguishing between phase cohesiveness and frequency synchronization. Note that 
phase cohesiveness can also be understood as a performance measure for synchroniza- 
tion and phase synchronization is simply the extreme case of phase cohesiveness with 
limt_j.oo 0{t) G ^(0). Indeed, if the magnitude r of the order parameter is understood 
as an average performance index for synchronization, then phase cohesiveness can 
be understood as a worst-case performance index. The following lemma relates the 
magnitude of the order parameter to a guaranteed level of phase cohesiveness. 

Lemma 2.1 (Phase cohesiveness and order parameter). Consider an array 
of n > 2 angles 9 ~ {6i, . . . ,9n) G T" and compute the magnitude of the order 
parameter r(9) = ^ \ X]j=i ^'^^ I- following statements hold: 

1) if 6 & A(7) for some 7 e [0,7r], then r{9) € [cos(7/2), 1]; and 

2) if r{9) G [0,1] and 9 e A(7r), then 9 G A(7) for some 7 e [2 arccos(r(6')), tt]. 
Proof. As customary, we abbreviate r{9) with r in what follows. The order 

parameter re^"^ is the centroid of all phasors e'*^ corresponding to the phases 9j when 
represented as points on the unit circle in C^. Hence, for 9 G A(7), 7 G [0,7r], it 
follows that r is contained in the convex hull of the arc of length 7, as illustrated 



in Figure 2.1 Let 7 G [0,7r] be fixed and let 6 G A(7). It follows from elementary 



Fig. 2.1. Schematic illustration of an arc of length -y S [0, tt], its convex hull (shaded), and the 
location • of the corresponding order parameter re^^ with minimum magnitude Tmin . 

geometric arguments that cos(7/2) — r,„in < r < r„iax = 1, which proves statement 
1). Conversely, if r is fixed and 9 G A(7r), then the centroid re"^ is always contained 
within the convex hull of the semi-circle A(7r) (centered at The smallest arc whose 
convex hull contains the centroid re"^ is the arc of length 7 — 2arccos(r) (centered 



at ■(/;), as illustrated in Figure 2.1 This proves statement 2). □ 

In the physics and dynamical systems community exponential synchronization is 
usually analyzed in relative coordinates. For instance, since the average frequency 



^ X]r=i ^i(^) = n Sr=i ^i ~ '^avg is coustaut, the Kuramoto model ( 1.1 1 is sometimes 



[49|, I38j analyzed with respect to a rotating frame in the coordinates £,i = 9i — Wavgi 
(mod 2tt), i E {1, . . . ,n}, corresponding to a deviation from the average angle. The 
existence of an exponentially stable one-dimensional (due to translational invariance) 
equilibrium manifold in ^-coordinates then implies local stability of phase-locked so- 
lutions and exponential synchronization. Alternatively, the translational invariance 



can be removed by formulating the Kuramoto model (1.1) in grounded coordinates 
Si = 9^ - 9nJoT i € {1, . . . ,n - 1} [5ni[3]. We refer to ppTLemma IV. 1] for a geomet- 
rically rigorous characterization of the grounded (5-coordinates and the relation of ex- 
ponential stability in (5-coordinates and exponential synchronization in ^-coordinates. 
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The following example of two oscillators illustrates the notion of phase cohesive- 
ness, applies graphical synchronization analysis techniques, and points out various 
important geometric subtleties occurring on the compact state space . 

Example 2.2 (Two oscillators). Consider n = 2 oscillators with W2 > loi. We 
restrict our attention to angles contained in an open half-circle: for angles 9i, 02 with 
1^2 — di\ < TT, we define the angular difference 02 — 0i to be the number in ] — 7r,7r[ 
with magnitude equal to the geodesic distance ~ ^i| and with positive sign iff the 
counter-clockwise path length from 0i to 02 on is smaller than the clockwise path 
length. With this definition the two-dimensional Kuramoto dynamics (0i,02) can be 
reduced to the scalar difference dynamics 02 — 0i- After scaling time as i H' i (w2 — wi) 
and introducing k = K/{uj2 — t^i) the difference dynamics are 

{O2 - 0i) = U02 - ^i) := 1 - Ksm{02 - 0i) . (2.1) 

at 

The scalar dynamics (|2.1[) can be analyzed graphically by plotting the vector field 



/k(^2 — ^1) over the difference variable 02 — 0i, as in Figure 2.2(a) Figure 2.2(a) dis- 
plays a saddle- node bifurcation at k = 1. For k < 1 no equilibrium of (2.1 ) exists, and 
for K > 1 an asymptotically stable equilibrium 6'stabie = arcsin(«;~^) e ]0, 7r/2[ together 
with a saddle point ^saddle = arcsin(K^^) e ]7r/2,7r[ exists. For 0{O) e A(|^saddic|) all 
trajectories converge exponentially to ^stable, that is, the oscillators synchronize expo- 
nentially. Additionally, the oscillators are phase cohesive iff 6{0) S A(|6'saddic|), where 
all trajectories remain bounded. For 0{O) ^ A(|6'saddio|) the difference 6*2 (i) — 0i{t) 
will increase beyond tt, and by definition will change its sign since the oscillators 
change orientation. Ultimately, 02{t) — 0i{t) converges to the equilibrium ^stable in 
the branch where 02 — 0i < 0. In the configuration space this implies that the dis- 
tance \02{t) — 0i{t)\ increases to its maximum value tt and shrinks again, that is, the 
oscillators are not phase cohesive and revolve once around the circle before converging 
to the equilibrium manifold. Since sin(0stabiG) = sin(6'saddic) — k^^, strongly coupled 
oscillators with k ^ 1 practically achieve phase synchronization from every initial 
condition in an open semi-circle. In the critical 1, the saddle point at 7r/2 

is globally attractive but not stable: for ^2(0) — 0i(O) =7r/2 + e (withe>0 sufficiently 
small), the oscillators are not phase cohesive and revolve around the circle before con- 



verging to the saddle equilibrium manifold in T^, as illustrated in Figure 2.2(b) Thus, 
the saddle equilibrium manifold is both attractor and separatrix which corresponds to 
a double zero eigenvalue with two dimensional Jordan block in the linearized case. 

In conclusion, the simple but already rich 2-dimensional case shows that two os- 
cillators are phase cohesive and synchronize if and only ii K > iiTcriticiai — uJ2~^i, and 
the ratio = Kcntidai/ K < 1 determines the ultimate phase cohesiveness as well 
as the set of admissible initial conditions. In other words, practical phase synchro- 
nization is achieved for K ^ -K^criticiai , and phase cohesiveness occurs only for initial 
angles 6(0) E A(7), 7 — arcsin(i^criticiai/-?^) & ]'''/2, 7r[. This set of admissible initial 
conditions A (7) can be enlarged to an open semi-circle by increasing if/i^criticiai- Fi- 
nally, synchronization is lost in a saddle- node bifurcatiorj^at K = /-Ccriticiai- In Section 
|4]we will generalize all outcomes of this simple analysis to the case of n oscillators. □ 

3. A Reviev^r of Bounds for the Critical Coupling Strength. In case that 
all natural frequencies are identical, that is, coi = u for all i G {1. . . . , n}, a transfor- 
mation to a rotating frame leads to w = 0. In this case, the analysis of the Kuramoto 



^For Kuramoto models of dimension n > 3, this loss of synchrony via a saddle-node bifurcation 
is only the starting point of a series of bifurcation occurring if K is further decreased, see I36| . 
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(a) Vector field (|2?TJ for 6»2 - 6»i > 



(b) Trajectory e{t) for k = 1 



Fig. 2.2. Plot of the vector field \2.1\ for various values of k and a trajectory 9{t) £ T'^ for the 
critical case k = 1, where the dashed line is the equilibrium manifold and ■ and • correspond to 
9(0) and limt_>oo S{t). The non- smoothness of the vector field f{92 — di) at the boundaries {0, n} is 
an artifact of the non- smoothness of the geodesic distance on the state space 



model is particularly simple and almost global stability can be derived by var- 
ious methods. A sample of different analysis schemes (by far not complete) includes 
the contraction property [35] . quadratic Lyapunov functions |29j . linearization [5], or 
order parameter and potential function arguments |43) . 

In the following, we review various analysis methods and the resulting bounds on 
the critical coupling strength for the case of non-identical frequencies. 

3.1. The Infinite Dimensional Kuramoto Model. In the physics and dy- 



namical systems communities the Kuramoto model (1.1) is typically studied in the 
continuum limit as the number of oscillators tends to infinity and the natural frequen- 
cies obey an integrable distribution function 17 : K — >■ K>o. In this case, the Kuramoto 
model is rendered to a first order continuity equation or a second order Fokker-Planck 
equation when stochasticity is included. For a symmetric, continuous, and unimodal 
distribution g{Lo) centered above zero, Kuramoto showed in an insightful and ingen- 
uous analysis [331 [M] that the incoherent state (i.e., a uniform distribution of the os- 
cillators on the unit circle) supercritically bifurcates for the critical coupling strength 

2 

-f^critical = 7777 ■ (3-1) 

7rg(0) 

The bound ( |3.1[ ) for the on-set of synchronization has also been derived by other 
authors, see |441 for further references. In [22] Ermentrout considered symmetric 
distributions g{Ld) with bounded domain uj g [— Wmax, '^max], and studied the existence 
of phase-locked solutions. The condition for the coupling threshold i^criticai necessary 
for the existence of phase-locked solution reads in our notation as [221 Proposition 2] 



max 



ifcritical peK,P>l 



\ C ^/p^^g[u)du\ . (3.2) 



Ermentrout further showed that formula (3.2 1 yields -ftTcriticai > 2winax for symmetric 



distributions and -ftTcriticai > ^'^max/T^ whenever g is non-increasing in [0,Wmax]- Both 
of these bounds are tight for a bipolar (i.e., a bimodal double-delta) distribution and 
a uniform distribution [22, Corollary 2], 148, Sections 3 & 4]. Similar results for the 
bipolar distribution are also obtained in [2, and in [37^ the critical coupling for a 
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bimodal Lorentzian distribution is analyzed. For various other references analyzing 
the continuum limit of the Kuramoto model we refer the reader to [44[ [5] . 

3.2. Necessary or Sufficient Bounds in the Finite Dimensional Ku- 
ramoto Model. In the finite dimensional case, we assume that the natural frequen- 
cies are supported on a compact interval G [wmaxi Wmin] C K for alH € {1, . . . , n}. 
This assumption can be made without loss of generality since the critical coupling 
^critical IS not finite for unbounded natural frequencies Wi gHl Theorem 1] . In [TH [55] 
a necessary condition for the existence of synchronized solutions states the critical 
coupling in terms of the width of the interval [wmax, '^min] as 

^> 2{n-l) ■ (^-^^ 

Obviously, in the limit as n — oo, this bound reduces (i^Jmax — <^min)/2, the simple 
bound derived in the introduction of this paper. A looser but still insightful necessary 
condition is K > 2a, where a is the variance of the uji [35], [JHl Corollary 2]. For 



bipolar distributions uji £ {i^min, Wmax}, necessary explicit conditions similar to (3.31 
can be derived for non-complete and highly symmetric coupling topologies [8] . 

Besides the necessary conditions, various bounds sufficient for synchronization 
have been derived including estimates of the region of attraction. Typically, these 
sufficient bounds are derived via incremental stability arguments and are of the form 

identical = ||'^C^||p-/K 7), (3.4) 

where is the p-norm and y is a matrix (of yet unspecified row dimension) mea- 
suring the non-uniformity among the cj^. For instance, V = In — (l/'^)lnxn gives 
the deviation from the average natural frequency, Vuj — uj — Wavglnxi- Finally, the 
function / : N x [0, 7r/2[— [1, oo[ captures the dependence of i^criticai on the number of 
oscillators n and the scalar 7 determining a bound on the admissible pairwise phase 
differences, which is, for instance, of the form ||(. . . , 9i{t) — 9j{t), ■ ■ ■)\\p < 7- 



Two-norm bounds, i.e., p = 2 in condition ( |3.4[ ), have been derived using quadratic 
Lyapunov functions in [Tli proof of Theorem 4.2] and [201 Theorem V.9], where 
the matrix V € K"'("-i)/2xn jg tj^g incidence matrix such that Vu; is the vector 
of n(7i — l)/2 pairwise differences LUi —ujj- A sinusoidal Lyapunov function [241 Propo- 
sition 1] leads to a two-norm bound with V — In — (l/'T')lnxn- Similar two-norm 
bounds have been obtained by contraction mapping |29[ Theorem 2] and by contrac- 
tion analysis (THl Theorem 3.8], where V G jjn-ixn orthonormal projector on 
the subspace orthogonal to l„xi- For all cited references the region of attraction is 
given by the n(n — l)/2 initial phase differences in two-norm or cx)-norm balls satis- 
fying ||T^0(O)|j2 00 Unfortunately, none of these bounds scales independently of 
n since ||l^aj||2 is a sum of at least n — 1 terms in all cited references and f{n,^) in 



condition (3.4) is either an increasing [29] or a constant function of n [20 | [T4| , [16 ] (24 



A scaling of condition (3.4) independently of n has been achieved only when 
considering the width Wmax — <^min = IK- ■ ■ it^i ~ 'j-'ji • ■ ■ )lloo' ^^^^ is, for Vuj being the 



vector of all n{n — l)/2 pairwise frequency differences and p = 00 in condition (3.4). 
A quadratic Lyapunov function leads to /(n, 7) — n/[2 sin(7)) [HI proof of Theorem 
4.1], a contraction argument leads to /(n, 7) — n/{{n — 2) sin(7)) [42] Lemma 9], and 
a geometric argument leads to the scale-free bound 7(7) = 1/(2 sin(7/2) cos(7)) [111 
proof of Proposition 1]. In [521 Theorem 3.3] and in our earlier work [201 Theorem V.3], 
the simple and scale-free bound /(7) = 1/ sin(7) has been derived by analyticity and 
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contraction arguments. In our notation, the region of attraction for synchronization 
is in ah cited references [HIH UHl [Ml HO] given as 61(0) G A(7) for 7 e [0,7r/2[. 

3.3. Implicit and Exact Bounds in the Finite Dimensional Kuramoto 
Model. Three recent articles [IHIEHIISI independently derived a set of implicit consis- 
tency equations for the exact critical coupling strength -ftTcriticai for which phase-locked 
solutions exist. Verwoerd and Mason provided the following implicit formulae to com- 
pute ^critical Theorem 3]: 

En , 
^i-{n,/u*Y, (3.5) 

where ilj = cj^ — ^ X^jLi ^ Ill^lloo ' ^ ll^llool unique solution to 

2^]" yr^oWF^^" i/v/r^oWF- (3.6) 

^ — ^z— 1 ^ ^l—l 

Verwoerd and Mason also extended their results to bipartite graphs [SD] but did not 



carry out a stability analysis. The formulae (3.5l-(3.6) can be reduced exactly to 
the implicit self-consistency equation derived by Mirollo and Strogatz in [38' and by 
Aeyels and Rogge in where additionally a local stability analysis is carried out. 
The stability analysis [351 13] in the n-dimensional case shows the same sadle-node 



bifurcation as the two-dimensional Example 2.2 for K < -fCcriticai there exist no 
phase-locked solutions, for K > iCcriticai there exist stable phase-locked solutions, and 
for K = i^criticai the Jacobians of phase-locked solutions equilibria have a double zero 
eigenvalue with two-dimensional Jordan block, as illustrated in Example [2. 2[ 

In conclusion, in the finite dimensional case various necessary or sufficient explicit 
bounds on the coupling strength -fCcriticai are known as well as implicit formulas to 
compute -fi'criticai which is provably a threshold for local stability. 

3.4. The Critical Coupling Strength for Second-Order Kuramoto Os- 
cillators. For rn — n, Di = 1, and Mi = M > the multi-rate Kuramoto model 



(1.3) simplifies to a second-order system of coupled oscillators with uniform inertia 
and unit damping. Such homogenous second-order Kuramoto models have received 
some attention in the recent literature [T3[|i^ liSlI^ l?7l [Tl [^ . 

In [13j two sufficient synchronization conditions are derived via second-order 



Gronwall's inequalities resulting in a bound of the form (3.4) with p = oo together 
with conditions on sufficiently small inertia or sufficiently large inertia [131 Theorems 
5.1 and 5.2]. In [T31 Theorem 4.1 and 4.2] phase synchronization was also found to 
depend on the inertia, whereas phase synchronization was found to be independent 
of the inertia in the corresponding continuum limit model [21 [1] . References [45l |46] 
observe a discontinuous first-order phase transition (where the incoherent state looses 
its stability) , which is independent of the distribution of the natural frequencies when 
the inertia M is sufficiently large. This result is also conffimed in [31[T]. In [25| a 
second-order Kuramoto model with time delays is analyzed, and a correlation between 
the inertia and the asymptotic synchronization frequency and asymptotic magnitude 
of the order parameter magnitude is observed. In [571 Ell] it is reported that inertia 
suppress synchronization for an externally driven or noisy second-order Kuramoto 
model, and [H [I] explicitly show that the critical coupling i^criticai increases with the 
inertia M for a Lorentzian or a bi-polar distribution of the natural frequencies. 

The cited results [T31 B51 \T7[ [H [2] on the inertial effects on synchronization 
appear conflicting. Possible reasons for this controversy include that the cited articles 
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consider slightly different scenarios (time delays, noise, external forcing), the cited re- 
sults are only sufficient, the analyses are based on the infinite-dimensional continuum- 



limit approximation of the finite-dimensional model ( 1.3 1, and some results stem from 



insightful but partially incomplete numerical observations and physical intuition. 

4. Necessary and Sufficient Conditions on the Critical Coupling. From 
the point of analyzing or designing a sufficiently strong coupling in the Kuramoto- 



type applications [331 13 |M1 HOI SZj , the exact formulae (3.5|-(3.6| to compute the 
critical coupling have three drawbacks. First, they are implicit and thus not suited 
for performance or robustness estimates in case of additional coupling strength, e.g., 
which level of ultimate phase cohesiveness or which magnitude of the order parameter 
can be achieved for K = c ■ i^criticai with a certain c > 1. Second, the corresponding 
region of attraction of a phase-locked equilibrium for a given K > -/^critical is unknown. 
Third and finally, the particular natural frequencies uJi (or their distributions) are 
typically time- varying, uncertain, or even unknown in the applications |44l [21 1431 1511 
\7U[ 137] . In this case, the exact ii'criticai needs to be dynamically estimated and re- 
computed over time, or a conservatively strong coupling K ^ i^criticai has to be chosen. 

The following theorem states an explicit bound on the coupling strength together 
with performance estimates, convergence rates, and a guaranteed semi-global region of 
attraction for synchronization. Besides improving all other bounds known to the au- 
thors, our bound is tight and thus necessary and sufficient when considering arbitrary 
distributions of the natural frequencies supported on a compact interval. 

Theorem 4.1. (Explicit, necessary, and sufficient synchronization con- 
dition) Consider the Kuramoto model ( |1.1[ ) with natural frequencies (wi, . . . ,w„) and 
coupling strength K . The following three statements are equivalent: 

(i) the coupling strength K is larger than the maximum non-uniformity among 
the natural frequencies, i.e., 

K > -/^critical — '^max ^min i 

(4.1) 

(ii) there exists an arc length 7max G ]'''/2,7r] such that the Kuramoto model ( |l.l[ ) 
synchronizes exponentially for all possible distributions of the natural frequen- 
cies supported on [wmin, '^max] and for all initial phases 9{Q) € A(7max),' o.'nd 



(Hi) there exists an arc length 7min € [0,7r/2[ such that the Kuramoto model (1.1 1 
has a locally exponentially stable synchronized trajectory in A(7inin) for all 
possible distributions of the natural frequencies supported on [wmin, i^max] • 
If the three equivalent cases (i), (ii), and (Hi) hold, then the ratio -R'criticai/^ dnd 
the arc lengths 7niin S [0; '''/2[ and 7niax S ]'''/2, tt] are related uniquely via sin(7uiin) = 
sin(7max) = -f^criticai/^j o,nd the following statements hold: 

1) phase cohesiveness: the set A(7) is positively invariant for every 7 S 
[7inin , 7max] 7 a'^'^ cach trajectory starting in A(7max) approaches asymptoti- 
cally A(7min); 

2) order parameter: the asymptotic value of the magnitude of the order pa- 
rameter denoted by rrx> — linit_>oo ^| X]j=i e'^^^*''|js bounded as 



1 > r„o > cos 



'1+ v/1- (^critical W 



3) frequency synchronization: the asymptotic synchronization frequency is 
the average frequency Wavg = ^X]"=i'^i' and, given phase cohesiveness in 
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A (7) for some fixed 7 < n/2, the exponential synchronization rate is no 
worse than Afg = X 005(7) ; ^'^'^ 
4) phase synchronization: if uji = s G R for all i G {!,..., n}, then for 
every 6{0) e ^(7), 7 € [0,7r[, the phases synchronize exponentially to the 
average phase 9avg{t) '■— ^ X]r=i ^C*^) + s • i (mod 2n) and the exponential 
synchronization rate is no worse than Aps = _R'sinc(7). 

To compare the bound (4.1) to the bounds presented in Section [sj we note from 
the proof of Theorem 4.1 that our bound (4.1 1 can be equivalently stated as if > 
(wmax ^ i^min)/ sin(7) and thus improves the sufficient bounds [HI [121 IHl IIHI HI] • 
In the simple case n — 2 analyzed in Example 2.2 the bound (4.1 1 is obviously exact 
and also equals the necessary bound (3.3). Furthermore, Theorem 4.1 fully generalizes 
the observations in Example |2.2| to the n-dimensional case. In the infinite-dimensional 
case the bound (4.1 1 is tight with respect to the necessary bound for a bipolar distri- 
bution oji € {wminji^max} derived in [^IHllllH]- Note that condition (4.1 1 guarantees 
synchronization for arbitrary distributions of Wi supported in [i^min , Wmax] , which can 
possibly be uncertain, time- varying (addressed in detail in Subsection 4.2 1, or even 
unknown. Additionally, Theorem |4.1| also guarantees a larger region of attraction 
61(0) G A(7,„ax) for synchronization than [lllEailSlliailHlEllEQiEe]. 



Besides the necessary and sufficient bound (4.1), Theorem 4.1 gives guaranteed 



exponential convergence rates for frequency and phase synchronization, and it estab- 
lishes a practical stability result in the sense that the multiplicative gap if critical /if in 
the bound (4.1) determines the admissible initial and the guaranteed ultimate phase 
cohesiveness as well as the guaranteed asymptotic magnitude r of the order parameter. 
In view of this result, the convergence properties of the Kuramoto model ( 1.1 ) are best 
described by the control-theoretical terminology "practical phase synchronization." 

The proof of Theorem 4.1 relies on a contraction argument in combination with 
a consensus analysis to show that (i) implies (ii) and thus also 1) - 4) for all natural 
frequencies supported on [wmim '^max]- In order to prove the implication (ii) =4> 
(i), we show that the bound (4.1) is tight: if (i) is not satisfied, then exponential 
synchronization cannot occur for a bipolar distribution of the natural frequencies. 
Finally, the equivalence (i), (ii) <^ (iii) follows from the definition of exponential 
synchronization and by basic arguments from ordinary differential equations 

Proof. Sufficiency (i) (ii): We start by proving the positive invariance 

of A(7), that is, phase cohesiveness in A(7) for some 7 g [0, tt]. Recall the geodesic 
distance on the torus and define the non-smooth function F : T" [0, tt]. 



Vitp) = maxflV-i - I i, j S {1, • ■ • 

The arc containing all initial phases has two boundary points: a counterclockwise 
maximum and a counterclockwise minimum. If we let /max(V-') (respectively imin(V')) 
denote the set indices of the angles ipi, . . . ,ipn that are equal to the counterclockwise 
maximum (respectively the counterclockwise minimum), then we may write 



y(V') = \-ipm' - i^l'l for all m G imax('0) and £' G /min(V')- 

By assumption, the angles 0i{t) belong to the set A(7) at time i = 0. We aim to 
show that they remain so for all subsequent times t > Q. Note that 9(t) € A(7) if and 
only if V{9{t)) < 7 < tt. Therefore, A(7) is positively invariant if and only if V{9{t)) 
does not increase at any time t such that V{9{t)) — 7. The upper Dini derivative of 
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V{9{t)) along the dynamical system (1.1) is given by (35] Lemma 2.2] 



where m € /max(^(i)) and £ G /min(^^(i)) are indices with the properties that 0m{t) = 
max{0,„/(t) I m' e /maxCfl*))} and 9i{t) = niin{^f.(i) | f e /min(6'(0)}- Written out 
in components D^V{9{t)) takes the form 

K " 

= _ - - ^(sin(0,„(t) - 0,{t)) + sm{e,{t) - ee{t))) . 

i=l 

Note that the index i in the upper sum can be evaluated for i G {l,...,rt}, and 
for i — m and i — i one of the two sinusoidal terms is zero and the other one 
achieves its maximum value in A(7). In the following we apply classic trigonometric 
arguments from the Kuramoto literature [Ml |42l [19]. The trigonometric identity 
sin(a::) + sin(y) = 2 sin (^±2^) cos leads to 

D^V{0{t)) = c.„, -.,-"^±(2 



^ .^,^t)^_OJt)-m^ (4.2) 



The equality V{6{t)) = 7 implies that, measuring distances counterclockwise and 
modulo additional terms equal to multiples of 2tt, we have 0m,{t) — 9i(t) — < 
d,n(t) ~ 9iit)<-f, and 0<9^{t) - 9i{t)<j. Therefore, D+V{9{t)) simplifies to 



D+Vi9{t)) < uj„, -uJi ^ (2 sin(^) cos 



1=1 

Reversing the identity from above as 2sin(a;) cos(?/) = sin(x — y) + sin(a; + y) yields 

K " 

D+V{9{t)) < w„ - W£ sin(7) = - w<> - K sin(7) . 

n ^ — ' 

2 = 1 

It follows that the length of the arc formed by the angles is non- increasing in A (7) if 
for any pair {to, f\ it holds that ii'sin(7) > ujm — ^jJi, which is true if and only if 

Ksin{-l)>K critical-, (4.3) 



where ii'criticai is as stated in equation ( [4.1| . For 76 [0, tt] the left-hand side of ( [4.3[ ) is 
a concave function of 7 that achieves its maximum at 7* = 7r/2. Therefore, there exists 
an open set of arc lengths 76 [0, tt] satisfying equation ([4.31) if and only if equation 



(4.3) is true with the strict equality sign at 7* = 7r/2, which corresponds to equation 



(4.1) in the statement of Theorem 4. 1 Additionally, if these two equivalent statements 



are true, then there exists a unique 7niin G [0,7i'/2[ and a 7niax G ]7''/2, tt] that satisfy 



equation (4.3) with the equality sign, namely sin(7min) — sin(7max) — -K'criticai/^- 
For every 7 € [7minj7max] it follows that the arc-length V{9{t)) is non-increasing, and 
it is strictly decreasing for 7 e ]7miii, 7max[- Among other things, this shows that 
statement (i) implies statement 1). 
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The frequency dynamics of the Kuramoto model ( 1.1 1 can be obtained by differ 



entiating the Kuramoto model (1.1 1 as 

d , 



dt 



Ell, . . 

_^^a,,{t)ie,~e,), (4.4) 



where aij{t) — (K/n) cos{9i{t) — Oj{t)). In the case that K > ii'criticai, we just proved 
that for every 9(0) G A(7max) and for all 7 G ]7min,7max] there exists a finite time 
T > such that 9{t) e A(7) for all t > T, and consequently, the terms aij{t) are 



strictly positive for all t > T. Notice also that system (4.4) evolves on the tangent 
space of T", that is, the Euclidean space K". Now fix 7 G ]7min,7r/2[ and let T > 
such that aij{t) > for all t>T, and note that the frequency dynamics (4.4) can be 
analyzed as the linear time-varying consensus system 

where L(t) = diag(^J_^j aij(t)) — A{t)) is a symmetric, fully populated, and time- 
varying Laplacian matrix corresponding to the graph induced by A{t). For each time 
instant t > T, the weights aij{t) are strictly positive, bounded, and continuous 
functions of time. Consequently, for each t > T the graph corresponding to L(t) 
is always complete and connected. Thus, for each t > the unique eigenvector 
corresponding to the zero eigenvalue is l„xi and l^^i^^ ^ 0- It follows that 
12^=1 ^i(t) — — "■'^avg is a conserved quantity. Consider the disagreement 

vector 5 = 9 — Wavglnxii as an error coordinate satisfying l^xi*^ = that is, S lives 
in the disagreement eigenspace of dimension n— \ with normal vector l„xi. Since Wavg 



is constant and ker(L(t)) = span(l„xi), the dynamics (4.4) read in (5-coordinates as 

±5 = -L{t)5. (4.5) 
Consider the disagreement function 5 !—>■ ||(5|p = (5^(5 and its derivative along the 



disagreement dynamics (4.5) which is ^ \\5\\^ = —2o^L{t)5. By the Courant-Fischer 



Theorem, the time derivative of the disagreement function can be upper-bounded 
(point-wise in time) by the second-smallest eigenvalue of the Laplacian L{t), i.e., 
the algebraic connectivity \2{L{t)), as ^ PIP < -2X2{L{t))\\8\\'^ . The algebraic 
connectivity \2{L{t)) can be lower-bounded as \2{L{t)) > ii'minj jgj-j^ „}{cos(0i — 
9j)\9 e A(7)} > K cos{j) = Afs . Thus, the derivative of the disagreement function 
is bounded as ^ \\S\\ < — 2Afs|j(5|p. The Bellman-Gronwall Lemma [321 Lemma A.l] 
yields that the disagreement vector 5{t) satisfies ||(5(t)|| < ||(5(0)||e~'*'f=' for all t > T. 
This proves statement 3) and concludes the proof of the sufficiency (i) (ii). 



Necessity (ii) (i): To show that the critical coupling in condition (4.11 

is also necessary for synchronization, it suffices to construct a counter example for 
which K < i^criticai and the oscillators do not achieve exponential synchronization 
even though all uji G [wmim'^max] and 9{0) G A(7) for every 7 G ]7r/2,7r]. A basic 
instability mechanism under which synchronization breaks down is caused by a bipolar 



distribution of the natural frequencies, as shown in Example 2.2 

Let the index set {1, . . . , n} be partitioned by the two non-empty sets Ii and l2- 
Let Ui = Wmin foi' « G Ii and uji = w^ax for i G I2, and assume that at some time 
t > it holds that 9i{t) = —7/2 for i G Ii and 6i{t) = +7/2 for i E I2 and for some 
7 G [0, 7r[. By construction, at time t all oscillators are contained in an arc of length 
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7 G [0,7r[. Assume now that /-C < ^'critical and the oscillators synchronize. Consider 
the evolution of the arc length V{9{t)) given as in (|4.2| by 



X cos 



ieii 

-,{t)-e,{t) 6,(t)-ei{t) 



2 2 



~ - |I (2 sin(^^^^^^^^^-^^ ^p„j^ ^m(t)-g»(0 ^ g.ffl-g^(t) ^^ 

where the summation is split according to the partition of {1, . . . , n} into Ii and 
By construction, we have that I ^Ii, m £ I2, oJi = Wmin, = ^max, Qi{t) — (^iit) = 
--ill for i e Ii, and B^ify = d„,(t) = +7/2 for i e la- Thus, DW(B{i)) simplifies to 

i?+V^(^^(t)) =C^„,ax-C^min-^ 51 S^^^ ( ^ ) ( ^ ) ) - ^ ^ (2 sin ( |) COS ^ ^ 



Again, we reverse the trigonometric identity via 2 sin(a:) cos(?/) = sin(a;— 2/)+sin(x+2/), 
unite both sums, and arrive at 

D^V(B{t)) = cj^,^ - cj^in - K sin(7) . (4.6) 

Clearly, for K < i^Tcriticai the arc length V{9{t)) = 7 is increasing for any arbitrary 
7 € [0, tt]. Thus, the phases are not bounded in A(7). This contradicts the assumption 
that the oscillators synchronize for K < i^criticai from every initial condition 9(0) € 
A(7). Thus, ifcriticai providcs the exact threshold. For K = i^criticai, wc know from 
[381 13] that phase-locked equilibria have a zero eigenvalue with a two-dimensional 
Jacobian block, and thus synchronization cannot occur. This instability via a two- 
dimensional Jordan block is also visible in (4.6) since D'^V(9{t)) is increasing for 



9{t) e A(7), 7 S ]7r/2, vr] until all oscillators change orientation, just as in Example 



2.2 This concludes the proof of the necessity (ii) (i). 

Sufficiency (i),(ii) =^ (iii)- Assume that (i) and (ii) hold and exponen- 



tial synchronization occurs. When formulating the Kuramoto model (1.1) in a ro- 
tating frame with frequency ujuvg, statement 3) implies exponential convergence of 
the frequencies 9i{t) to zero. Hence, for all 9(0) G A(7inax) every phase 9i{t) con- 
verges exponentially to a constant limit phase given by ^i^sync = linif-i.00 (^) = 
9i{0) + 9i{T), which corresponds to an equilibrium of the Kuramoto model ( |1.1[ ) 
formulated in a rotating frame. Furthermore, statement 1) implies that these equi- 
libria (6'i,sync, • ■ • , ^n,sync) are Contained in A(7min). Since the Kuramoto model ( |1.1[) 
features only a finite number of fixed points (modulo translational invariance) [6l 
Lemma 1.1 and Theorem 4.1], almost every trajectory with 0(0) € A(7i„ax) converges 
exponentially to a one-dimensional stable equilibrium manifold. Hence, if condition (i) 
holds, there exists a locally exponentially stable synchronized solution 9{t) G A(7i„in)- 
Necessity (iii) (i),(ii): Conversely, assume that condition (i) does not 

hold, that is, K < ^critical — i^max ~ i^min. Wc provc the necessity of (i) again by 
invoking a bipolar distribution of the natural frequencies. In this case, it is known that 
for K — i^criticai = "^max ~ '^min there cxists a unique equilibrium (in a rotating frame 
with frequency Wavg), and for K < Kcnticai there exists no equilibrium |48l Section 4]. 
In the latter case, synchronization cannot occur. In the former case, the equilibrium 
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configuration corresponds to the phases arranged in two clusters (sorted according 
to the bipolar distribution) which are exactly 7r/2 apart [151 Section 4]. Finally, 



note that such an equilibrium configuration is unstable, as shown by identity (4.6) 



We remark that the same conclusions can alternatively be drawn from the implicit 
equations (3.5|-(3.6| for the critical coupling. This proves the necessity (iii) ^ (i)'(ii)- 

By statement 1), the oscillators are ultimately phase cohesive in A(7niin). It fol- 
lows from Lemma [2.1| that the asymptotic magnitude r of the order parameter satisfies 
1 > r > cos(7„iin/2). The trigonometric identity cos(7i„in/2) = ^^(1 + cos(7inm))/2 
together with a Pythagorean identity yields the bound in statement 2). 

In case that all natural frequencies are identical, that is, uji — s for all i G 
statement 1) implies that 7,nin — and 7max t '^^ In short, the phases 
synchronize for every 6{Q) G A(7r). The coordinate transformation 9 i-^ 9 + st yields 
the dynamics 9i = — bij{t){9i — Oj), where bij{t) = {K/n)sinc{9i{t) — 9j{t)) is 

strictly positive for all i > due to the positive invariance statement 1). Statement 
4) can then be proved along the lines of statement 3). □ 

4.1. Statistical studies. Theorem |4 . 1 1 places a hard bound on the critical cou- 
pling strength i^criticai for all distributions of Wi supported on the compact interval 



This set of admissible distributions includes the worst-case bipolar distri- 

For a particular distribution g{u;) supported 



bution used in the proof of Theorem 
on 



4.1 



the bound (4.11 is only sufficient and possibly a factor 2 larger than 



the necessary bound (3.3 1. The exact critical coupling for g{Lo) lies somewhere in 



between and can be obtained by solving the implicit equations (|3.5|)-(3.6) 



The following example illustrates the average case for natural frequencies sampled 
from a uniform distribution g{Ld) = 1/2 supported for uj E [—1, 1]. Figure 4.1 reports 
numerical findings on the critical coupling strength for n e [2, 300] oscillators in a 
semi-log plot, where the coupling gains for each n are averaged over 1000 simulations. 




Fig. 4.1. Analysis of the necessary bound | |3.3| 
sufficient explicit bound (|4.1[l (u) 



(O), the exact bound ||3.5|-([3^ (0), and the 



First, note that the three displayed bounds are equivalent for n = 2 oscillators. 



As the number of oscillators increases, the sufficient bound (4.1 1 clearly converges to 
~ '^min — 2, the width of the distribution (/(w), and the necessary bound (3.31 



accordingly to half of the width. The exact bound (3.5)-(3.6) quickly converges to 



16 



F. Dorfier and F. BuUo 



4(wmax ~ '^min)/(27r) — A/iT in agreement with the results (3.1) and (3.2) predicted in 
the case of a continuum of oscillators. It can be observed that the exact bound (3.5 )- 



(3.6) is closer to the sufficient and tight bound (4.1 ) for a small number of oscillators, 

For large n, the 
In this case, the 



i.e., when there are few outliers increasing the width Wmax — i^mii 
sample size of uji increases and thus also the number of outliers 



exact bound (3.5)-(3.6) is closer to the necessary bound (3.3) 



4.2. Extension to time-varying natural frequencies. One motivation to 



prefer the explicit and tight bound (4.1 1 over the implicit and exact bound (3.5)-(3.6l 
are time- varying natural frequencies uji{t) bounded in [wmax: "^min]- We distinguish 
the two cases of switching and slowly and smoothly time-varying natural frequencies 



and note that Theorem 4.1 and its proof can be easily extended to these cases. 



4.2.1. Piece-wise constant uji{t). Consider a sequence of time instances {ifcjfegN 
such that to = and t^+i > tk for all fc g N. Assume that the natural frequencies 
U!i{t) are constant and bounded in [wmax, i^min] within each interval t G [tk,tk+i[- At 
time-point t^+i the natural frequencies may be discontinuous and switch. Note that 
the synchronization frequency and the corresponding phase-locked equilibria (on a 
rotating frame) will change with every switching instant. 

In this case, between any two switching instants, t G [tk,tk+i[, our analysis still 
holds and Theorem |4.1| can be applied without any modification. For all time t > 
and for all 6* e A (7), 7 e ]7inin, Tmax]: the arc length V(^9{t)) is strictly and uniformly 
decreasing for any switching sequence {tfcjfcgN, i-e, it is a so-called common Lyapunov 
function. As an outcome, the ultimate phase cohesiveness in A^-j^i^) will always 
be reached asymptotically despite the switching natural frequencies. Furthermore, if 
there exists a uniform dwell time e > such that t^+i —tk>e for all fc G N, then the 
derived synchronization rate Afs admits an estimate on \\in.f\tk+i II ^(^) ^'^avg(i)lnxi||j 
that is, how close the oscillators come to frequency synchronization within each in- 
terval [ti;,ifc+i[. Figure 4.2 illustrates all of these conclusion in a simulation. 




Fig. 4.2. Simulation of a network of n — 10 Kuramoto oscillators satisfying ./^/i^critical — 

1.1, 

where the natural frequencies aji(t) and LO„(t) (displayed in red dashed lines) are switching between 
constant values in [ojmin , '^max] = [0,1]. The simulation illustrates the phase cohesiveness of the 
angles 9{t) in A(7,jiiii), the exponential convergence of the freguencies 9(t) towards i.Liavg(t) between 
consecutive switching instances, as well as the monotonicity ofV{9(t)) in A{-y) for ^ S [7min, 7max]. 



In comparison, the analysis schemes [49l |3l |38] have to re-compute the exact 



implicit bound (3.5)-(3.6) after every switching instant, since they explicitly make use 



of the values of oJi and the corresponding equilibria. Obviously, the analysis schemes 
[iniinilSH] fail entirely in the case of time- varying frequencies analyzed in the following. 
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4.2.2. Slowly and smoothly varying Wiit). For smooth functions bounded 
Wmaxi^niin], the proof of phasB cohesiveness can be adapted without major modi- 



m 

fications 



However, the Kuramoto frequency dynamics (4.4 1 are rendered to 



— I 



E 



i=i 



(4.7) 



where aij(t) 



(K/n) cos {9 i{t) — Oj{t)) as before. The forced frequency dynamics 



(4.7) can be analyzed on the subspace orthogonal to 1 
varying disagreement vector 6(t) = 0(t)— Wavg(i)l 



1 by considering the time- 
1 , as an error coordinate satisfying 



^nxi^i^) — 0- The frequency dynamics (4.7 1 read then in ^-coordinates as 



dt 



5 = n{t) - L{t)5, 



(4.8) 



where Cl{t) = uj{t) — a;avg(i)lnxi- On the subspace orthogonal to l„xi the dynamics 
(4.8) are exponentially stable for Q.{t) = 0, and a time-varying equilibrium frequency 
can be uniquely obtained as 5{t) = L\t)Vl{t), where is the Moore- Penrose inverse 
of L. In this case, the standard theory of slowly varying systems [Ml Chapter 9.6] can 
be applied for a slowly varying (l{t) satisfying ||n(t)||oo < e for e sufficiently small. 

In summary, ifeacha;i(<) is a smooth, bounded in [wmax, i^min], and the relative ac- 
celeration ||r2(i)||oo = ||w(t)— a)avg(i)lnxi||oo < £ is sufficiently small, then there exists 
T > and k = fc(e) > such that the frequencies satisfy ||<5(t) — {t)Vl{t)\\oo < k for 
all t>T. Moreover, if VL[t) — >• as t — >■ oo, then 5{t) -> V' {t)h{t) as t — >■ oo. In par- 
ticular, e and k depend on the phase cohesiveness (5(t) G A(7), see [SUl Theorem 9.3] 
for details. Figure [473| illustrates these conclusions with a simulation of n = 10 oscil- 
lators. The authors of [MIIIH] come to a similar conclusion when analyzing the effects 
of time- varying frequencies via input-to-state stability arguments or in simulations. 





3 

^2.5 
^ 2 
Sl.5 
^ 1 
0.5 



t\s] 



10 



Fig. 4.3. Simulation of a network of n = 10 Kuramoto oscillators satisfying ^/-K^critical = 1-1; 
where the natural frequencies uji : K>o — ^ [i^min i i^max] = [0,1] are smooth, bounded, and distinct 
sinusoidal functions. Ultimately, each natural frequency u)i(t) converges to uji + sin(7rt) with uii £ 
[0, 1], and thus the relative acceleration Q{t) = u){t) — Li;avg(t)lnxl converges to zero. The simulation 
illustrates the phase cohesiveness of the angles 9{t) in A(7iiiin); iirid the ultimate boundedness of 
the frequency variations (disagreement vector) 5{t) = 8(t) — a;avg(i)lnxl ci,nd their convergence to 
zero. The simulation further confirms the monotonicity of V{0{t)) in A(7) for ■y S [7min, 7max]. 
Ultimately, V{9{t)) converges to a constant value (strictly below 7niinj o-s the frequencies converge. 



5. Synchronization of Multi-Rate Kuramoto Models 

extend the results in Theorem 



4.1 



In this section we 
to the multi-rate Kuramoto model (1.3 1. For the 
special case of second-order oscillators (m = n) with unit damping Di — 1 and uniform 
inertia Mi ^ M > 0, the literature US m UHl 123 H 12] on the inertial effects on 
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synchronization is controversial. Here we will rigorously prove that the inertial terms 
do not affect the location and local stability properties of equilibria of the multi-rate 
Kuramoto model (1.3). In particular, the necessary and sufficient synchronization 
conditions as well as the synchronization frequency are independent of theinertiae Mi] 
they rather depend on the terms Di mimicking viscous damping and time constants. 

5.1. A One-Parameter Family of Dynamical Systems and its Proper- 



ties. In this subsection we will link the multi-rate Kuramoto model ( 1.3 ) and the first- 



order Kuramoto model (1.1) through a parametrized system. Consider for ni,n2 > 
and A G [0, 1] the one-parameter family "Hx of dynamical systems combining dissipa- 
tive Hamiltonian and gradient-like dynamics together with external forcing as 

DiXi=Fi-\/iH{x), 





'In, ■ 




X2 


nx : 


M 




±3 



(1 - A)F2 
(1-A) 



(5.1) 



/„; 




AD^^ ■ 




'V2H{x) 







o" D2 


) 


V3-ff(x)_ 



where x = (xi, X2, X3) G Xi x X2 ^ K"^ = X is the state, and the sets Xi and X2 are 
smooth manifolds of dimensions ni and 712, respectively. The matrices Di g 
D2 E M"2xn2 g^j^fj ][f ^ ]]jn2xn2 g^j-g positive definite, are zero matrices of appropriate 
dimensiorj^ Fi e lR"i and F2 G M"^ are constant forcing terms, and H : A" — M is 
a smooth potential function with partial derivative \IiH{x) = dH{x)/dxi, gradient 
vector VH{x) = {dH(x)/dx)'^ G and Hessian matrix V^H{x) G M"""". 



The parameterized system (5.11 continuously interpolates, as a function of A G 
[0, 1], between gradient-like and mixed dissipative Hamiltonian/gradient-like dynam- 
ics. For A = 1, the system (5.11 reduces to gradient-like dynamics with forcing term 
F = [F^ , F2 , 0]"^ and time constant (or system metric) D = blkdiag(Z3i, 1)2, £'2^^-^) as 



-Hi : Di: = F - VH{x) . 



(5.2) 

For A = 0, the dynamics ( |5.1[ ) reduce to gradient-like dynamics for xi and dissipative 
Hamiltonian (or Newtonian) dynamics for {x2,X3) written as 



Dix, - Fi - ViH{x) , 



Ho : 



'In, ■ 




X2 









/„; 










) 


'V2H{x) 


M 




.^3. 








-/„2 







D2 




V^Hix) 



(5.3) 



It turns out that, independently of A G [0, 1], all parameterized systems of the form 
(5.1) have the same equilibria with the same local stability properties determined by 
potential function H{x). The following theorem summarizes these facts. 

Theorem 5.1 (Properties of Hx family). Consider for A G [0,1] the one- 
parameter family Hx of dynamical systems (5.1). The following statements hold in- 



dependent of the parameter A G [0, 1] and independent of the particular positive definite 
matrices Z?i, D2, M : 

1. Equilibria: For all A G [0, 1] the equilibria ofHx are given by the set £ = 
{x e X : VH{x) ^ F}; and 

2. Local stability: For any equilibrium x* € £ and for all A G [0,1], the 
inertia of the Jacobian of Hx is given by the inertia of —'S7^H{x*) and the 
corresponding center- eigenspace is given by the nullspace of'S7^H{x*). 



^We did not index tlic zero matrices according to their dimension to avoid notational clutter. 
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Statements 1) and 2) assert that normal hyperbolicity of the critical points of 
H{x) can be directly related to local exponential (set) stability for any Ae [0, 1]. This 
implies that all vector fields Ha, A e [0, 1], are locally topologically conjugate |39J near 
a hyperbolic equilibrium point e In particular, near x* Cz £, trajectories of 



the gradient vector field (5.2) can be continuously deformed to match trajectories of 
the Hamiltonian vector field (5.3) while preserving parameterization of time. This 
topological conjugacy holds also for hyperbolic equilibrium trajectories [171 Theorem 
6] considered in synchronization. The similarity between second-order Hamiltonian 
systems and the corresponding first-order gradient flows is well-known in mechanical 
control systems [511 [5^ , in dynamic optimization [H O US] , and in transient stability 
studies for power networks [12l [11] [15] , but we are not aware of any result as general as 
Theorem 5.1 In [IH[IIl[in], statements 1) and 2) are proved under the more stringent 
assumptions that "H^ has a finite number of isolated and hyperbolic equilibria. 

REMARK 5.2 (Extensions on Euclidean state spaces). If the dynamical sys- 
tem 7i\ is analyzed on the Euclidean space then it can be verified that 
the modified potential function H : i?"i+2n2 _^ jj{x) = -F[xi - ^2^2:2 + 
H{xi,X2, M^/^x^) is non-increasing along any forward-complete solution x : M>o — >■ 
^ni+2«2 and for all A e [0,1]. Furthermore, if the sublevel set = {x & X : 
H{x) < c} is compact, then every solution initiating in flc is bounded and forward- 
complete, and by the invariance principle (SO] Theorem 4.4] it converges to the set 
£ n Sic, independently of A S [0,1]. These statements can be refined under further 
structural assumptions on the potential function H{x), and various other minimizing 
properties can be deduced, see [11[51[1S]. Additionally, if H{x) constitutes an energy 
function, if all equilibria are hyperbolic, and if a one-parameter transversality condi- 



tion is satisfied, then the separatrices of system (5.1 1 can be characterized accurately 
[121 [m [15] . For zero forcing F = 0, these convergence statements also hold on the 
possibly non-Euclidean space X, and for non-zero forcing they hold locally on X. □ 



Proof. To prove statement 1), we reformulate the parameterized dynamics (5.1 1 as 



Xi 




"L>fi 







'Fi - ViH{x 


±2 




XD^^ 


-(l-A)/„, 




F2 - V2H{x 


_Mi3_ 




(1 - A)/„, 


D2 




-V^Hix) 















It follows from the Schur determinant formula 52] that det(WA) = dei{D^^){X+ (1 — 
A)^) is positive for all A G [0, 1]. Hence, W\ is nonsingular for all A G [0, 1], and the 
equilibria of (5.1) are given by by the set £ — {x € X : VH{x) = F}. To prove 
statement 2) we analyze the Jacobian of Hx at an equilibrium x* €z £ given by 






















XD-^ 


(A 


- l)M-i 









(1 - X)M-' 


AF 


-^D2M-^ 





Ini+n2 











V^Hix*). (5.4) 



Again, we obtain det(S'A) dct(i:ij^^) dct(i:)^^) det(Af"ii:>2M"^)(A+(l-A)2). Thus, 
is nonsingular for A S [0, 1], and the nuUspace of Jx{x*) is given by kerW^H{x*) 
(independently of A € [0, 1]). To show that the stability properties of the equilibrium 
X* & £ are independent of A e [0, 1], we prove that the inertia of the Jacobian J\{x*) 
depends only on S{x*) and not on A € [0, 1]. For the invariance of the inertia we 
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appeal to the main inertia theorem for positive semi-definite matrices [TOl Theorem 
5]. Note that J\{x*) and J\{x*Y' have the same eigenvalues. Let A = Jx{x*)'^ and 
P = S{x*)^ and consider the matrix Q defined via the Lyapunov equality as 



Q^-{AP + PA^) = P 


















P. 



Note that Q is positive semidefinite for A > 0, and for A 7^ the nuUspaces of Q and 
P coincide, i.e., kerQ ~ kerP. Hence, for A e ]0, 1] the assumptions of [101 Theorem 
5] are satisfied, and it follows that the non-zero inertia of A — Jx{x*)^ (restricted 
to image of A) corresponds to the non-zero inertia of P. Hence, the non-zero inertia 
of J\{x*) is independent of A G ]0, 1], and possible zero eigenvalues correspond to 
kerJA(a;*) — kerV^i?(a;*). To handle the case A = we invoke continuity arguments. 
Since the eigenvalues of J\{x*) are continuous functions of the matrix elements, the 
inertia of Jq{x*) is the same as the inertia of J\{x*) for A > sufficiently small. Since 
the inertia of J\{x*), X G ]0, 1], equals the inertia of P (which is independent of A), it 
follows that the inertia of J\{x*) equal the inertia of P for all A G [0, 1]. 

Finally, since blkdiag(/„j4.„2 , M) is positive definite, Sylvester's inertia theorem 
[TU] asserts that the inertia of P = blkdiag(/„j+„2 , M)(— V^i/(a;*)) equals the inertia 
of —W^H{x*). In conclusion, the inertia and the nuUspace of J\{x*) e qual the inertia 
of -V^i/(x*) and kerV^i? (x*). This completes the proof of Theorem 



5.1 



□ 



5.2. Equivalence of Local Synchronization Conditons. As a consequence 



of Theorem 5.1 we can link synchronization in the multi-rate Kuramoto model (1.31 



and in the regular Kuramoto model ( 1.1 ). Since Theorem 5.1 is valid only for equilib 



ria, we convert synchronization to stability of an equilibrium manifold by changing co- 
ordinates to a rotating frame. The explicit synchronization frequency Wgync G M of the 
multi-rate Kuramoto model ([1.3^ is obtained by summing over all equations (1.3 1 as 



(5.5) 



In the frequency-synchronized case when all 9i — and 6i — Wsync, equation (5.5) sim- 
plifies to J2i=i DiLUsync — X]r=i ^i- conclude that the synchronization frequency of 
the multi-rate Kuramoto model is given by Wgync — T^^^i^il T^l^i-^i- Accordingly, 
define the first-order multi-rate Kuramoto model by dropping the inertia term as 

K 



DiOi = I 



Ell 
1- 



sm{Qi - 0j) , i G {1, . . . ,71} , 



and the globally exponentially stable frequency dynamics as 



dt 



'^sync 



), «e{i,. 



A, 



(5.6) 



(5.7) 



where Mi, Di, Wj, and K take the same values as the corresponding parameters for 



the multi-rate Kuramoto model (1.3). It can be verified that the multi-rate Kuramoto 



model (1.3) and its first-order variant (5.6) have the same synchronization frequency. 
Finally, let uii = uii — DiUJsync and define the scaled Kuramoto model by 



K 



■ LUi 



j), iG{l,...,n}, 



(5.8) 
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and its associated scaled frequency dynamics by 



,to}. 



(5.9) 



The scaled model (5.8)-(5.9) corresponds to the dynamics (5.6)-(5.7) formulated in a 



rotating frame with frequency Wsync and after unitizing the time constants Di in (5.6) 



Notice that the multi-rate Kuramoto model (1.3), its first-order variant (5.61 



together with frequency dynamics ( 5.7 ) (formulated in a rotating frame with frequency 
Wsync), and the scaled Kuramoto model (|5.8|) together with scaled frequency dynamics 



(5.9) are instances of the parameterized system (5.1) with the forcing terms uji and 

" ^-YEL=iCos( 



the potential : T" x M" ^ M, H{e, 9) = \e'^ 0-^^ E"j=i cos(6l, -6*^) defined up to 
a constant value. In the sequel, we seek to apply Theorem |5.1| to these three models. 
For a rigorous reasoning, we define a two-parameter family of functions (j)r,s '■ 
T of the form 4'r,s{t) ^ r + s ■ t (mod 27r), where r e T and s e M. Consider 



for (ri. 



e A(7), 7 e [0,7r[ the composite function 



7,^ 



s{t), ■ ■ • ,0r„,s(i)) 



(5.10) 



mimicking synchronized trajectories of the three Kuramoto models (|1.3|, (5.6), and 



(5.^ 



We now have all ingredients to state the following result on synchronization. 
Theorem 5.3. (Synchronization Equivalence) Consider the multi-rate Ku- 



ramoto model (|1.3|), its first- order variant (|5.6|) 
with ui. 



LUi - DiUJs 



where ujg 



and the scaled Kuramoto model (5.81 
J2k=i'^k/ J2k=i ■ The following statements 



are equivalent for any 7 G [0,7r[, t >0, and any function ^*7,a;ay„c (*) defined in (5.10); 
(i) {^j^i^^^^^{t),(jJsync'^mxi) paramctrizcs a locally exponentially stable synchro- 



nized trajectory (9(t),6{t)) of the multi-rate Kuramoto model (1.3); 
^■y,u:sync W parametrizes a locally exponentially stable synchronized trajectory 
9{t) of the first- order multi-rate Kuramoto model (5.6); and 



(n) 

(Hi) $^,o(0 parametrizes a locally exponentially stable synchronized equilibrium 



trajectory 9{t) of the scaled Kuramoto model (5.8) 



If the equivalent statements (i), (ii), and (Hi) are true, then, locally near their respective 



synchronization manifolds, the multi-rate Kuramoto model (1.3 I, its first- order vari- 
ant (5.6) together with the frequency dynamics (5.7), and the scaled Kuramoto model 



(5.8) together with the scaled frequency dynamics (5.9) are topologically conjugate. 
For purely second-order Kuramoto oscillators (|1.3[) (with n = to). Theorem 5.1 



and Theorem |5.3| essentially state that the locations and stability properties of the 
foci of second-order Kuramoto oscillators (with damped oscillatory dynamics) are 
equivalent to those of the nodes of the scaled Kuramoto model (5.8) and the scaled 
frequency dynamics (|5.9[) (with overdamped dynamics), as illustrated in Figure 5.1 



Proof. By Definition, a synchronized trajectory of the multi-rate Kuramoto model 

^^(i),tJsynclmxi) for 7 e [0, 7r[ and i > 0. In 



(1.3) is of the form {9{t),9{t)) £ ($^ 



a rotating frame with frequency Wgync the multi-rate Kuramoto model (1.3) reads as 



K 



M9i -I- DiOi ^ uji > sin{9i ~ 9^) , i e {1, . . . , to} . 



_ K ^ 

Di9i = LOi } sin( 



(5.11) 



9j) , i G {to + 1, . . . , n} . 



Hence, an exponentially synchronized trajectory of (5.11) is an equilibrium solu- 
tion determined up to a translational invariance in §^ and satisfies {9(t),9(t)) € 
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Fig. 5.1. Phase space plot of a network of n = A second-order Kuramoto oscillators with 
n = m (left plot) and the corresponding first-order scaled Kuramoto oscillators | |5.8[ l together with 
the scaled frequency dynamics ^5.9\ (right plot). The natural frequencies uii, damping terms Di, 
and coupling strength K are such that uJsync = and K / K^,.\i\^^i = 1.1. From the same initial 
configuration 9(0) (denoted by U) both first and second-order oscillators converge exponentially to 
the same nearby phase-locked equilibria (denoted by *) as predicted by Theorems\5 . 1\ artd|5. 3[ 



($.y_o(i), Otoxi)- Hence, the phase-synchronized orhit ($.y_o(*), Omxi), understood as 
a geometric object in T" x K™, constitutes a one-dimensional equihbrium manifold of 
the multi-rate Kuramoto model (5.111. After factoring out the translational invari- 
ance of the angular variable 0, the exponentially-synchronized orbit (<&7,o(Oj Omxi) 
corresponds to an isolated equilibrium of (5.11) in the quotient space T" \ x M™. 



Since an isolated equilibrium of a smooth nonlinear system with bounded and Lips- 
chitz Jacobian is exponentially stable if and only if the Jacobian is a Hurwitz matrix 
[301 Theorem 4.15], the locally exponentially stable orbit ($7,o(i), Omxi) must be 
hyperbolic in the quotient space T" \ x M™. Therefore, the equilibrium trajec- 
tory ($^^o(i)i Omxi) is exponentially stable in T" x M™ if and only if the Jacobian of 
(5.11) evaluated along ($7,0(0:0^x1), has n-|-m — 1 stable eigenvalues and one zero 
eigenvalue corresponding to the translational invariance in 

By an analogous reasoning we reach the same conclusion for the first-order multi- 



rate Kuramoto model (5.6) (formulated in a rotating frame with frequency Wgync) 
and for the scaled Kuramoto model ( |5.8[ ): the exponentially-synchronized trajectory 
^^fl{t) e T" is exponentially stable if and only if the Jacobian of (5.8) evaluated 
along ^^fi{t) has n — 1 stable eigenvalues and one zero eigenvalue. Finally, recall 
that the multi-rate Kuramoto model (5.11), its first-order variant (5.6) together with 



frequency dynamics (5.7) (in a rotating frame), and the scaled Kuramoto model (5.8) 
together with scaled frequency dynamics ( |5.9[ ) are all instances of the parameterized 
system (5.1 ). Therefore, by Theorem |5.1[ the corresponding Jacobians have the same 
inertia and local exponential stability of one system implies local exponential stability 
of the other system. This concludes the proof of the equivalences (i) (ii) (iii). 
We now prove the final conjugacy statement. By the generalized Hartman- 
1 [171 Theorem 6], the trajectories of the three vector fields (5.11), 



Grobr 



the 



(5.6)-(|5.7[) (formulated in a rotating frame), and (5.8)-(5.9) are locally topologically 



conjugate to the flow generated by their respective linearized vector flelds (locally 
near {^j,o{t),Onixi)- Since the three vector fields ( |5.11[ ), ([5!6|-([5J1), and ([5^-([5!9| 
are hyperbolic with respect to {^j^oit), Omxi) and their respective Jacobians have the 
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same hyperbolic inertia (besides the common one-dimensional center eigenspace cor- 
responding to {^-yfi(t),Omxi)j the corresponding three linearized dynamics are topo- 
logically conjugate 1395 Theorem 7.1]. In summary, the trajectories generated by the 



three vector fields (5.11), (5.6)-([5T7| (formulated in a rotating frame), and (5.8)-(5.9l 



are locally topologically conjugate near the equilibrium manifold ($^.0(^)1 0„jxi)- D 

REMARK 5.4 (Alternative ways from first to second-order Kuramoto 
models) . Alternative methods to relate stability properties from the first-order Ku- 
ramoto model ( |1.1[ ) to the multi-rate model (1.31 include second-order Gronwall's 
inequalities |13j , strict Lyapunov functions for mechanical systems |311 132) , and sin- 
gular perturbation analysis [21] . It should be noted that the approaches |T3l |3ll [32] 
are limited to purely second-order systems, the second-order Gronwall inequality ap- 
proach |13j has been carried out only for uniform inertia = M and unit damping 
Di = 1, and the Lyapunov approach [311 132j is limited to potential-based Lyapunov 
functions and seems not extendable to our contraction-based Lyapunov function used 
in the proof of Theorem 4.1 Finally, the singular perturbation approach [2I| requires 
a sufficiently small inertia over damping ratio e = max^g^i^ ,„j.{Mi/Z)i}. 

As compared with these alternative methods. Theorem |5.3| applies to the multi- 



rate Kuramoto model (1.3) with mixed first and second-order dynamics, for all values 
of Mi > and Di > 0, and without additional assumptions. Finally, it is instructive 



to note that the first-order multi-rate Kuramoto dynamics (5.6) and the frequency dy- 



namics (5.7 1 (in the time-scale t/e) correspond to the reduced slow system and the fast 
boundary layer model in the singular perturbation approach |2H Theorem IV. 2]. □ 

5.3. Synchronization in the Multi-Rate Kuramoto Model. Theorems |5.1| 
and 5.3 together with Theorem 4.1 on the first-order Kuramoto model ( 1.1 1 allow us to 



state our final conditions on synchronization in the multi-rate Kuramoto model (1.3). 

Theorem 5.5. (Exponential Synchronization in the Multi-Rate Ku- 
ramoto Model) Consider the multi-rate Kuramoto model (1.3) and let Coi = Ui — 
DiUisync, where Wgync = X]fc=i ^k/ X]fc=i ^k- Then the following statements hold: 
1) Exponential synchronization: The following two statements are equivalent: 
(i) the coupling strength K is larger than the maximum non-uniformity among 

the scaled natural frequencies, i.e., K > -fsTcriticai — '^max — ^^mini (ind 
(ii) there exists an arc length 7min G [0,7r/2[, such that the multi-rate Kuramoto 



model (1.3) has a locally exponentially stable synchronized solution with syn- 
chronization frequency Wgync 0,'rid phase cohesive in A(7min) for all n > 2, for 
all m e [0,n], for all inertiae Mj > 0, j € {1, . . . and for all possible 

A > anduii £ R satisfying uJi = w — A^^sync e [wmax, Wmin], i G {1, . . . ,n}. 
Moreover, in either of the two equivalent cases (i) and (ii), the ratio i^criticai/-?^ o,nd 
the arc length 7min G [0i''''/2[ are related uniquely via Kcriticai/ K ~ sin(7,nin)- 
2) Phase synchronization: The following two statements are equivalent: 
(Hi) there exists a constant s G M such that uJi — DiS for all i £ {\, . . . ,n} ; and 
(iv) there exists an almost globally exponentially stable phase-synchronized solu- 



tion with constant synchronization frequency oj. 



sync 



Moreover, in either of the two equivalent cases (Hi) and (iv), the constant s and the 
synchronization frequency Wgync o,re related uniquely via s = Wsynci o.nd the the asymp- 
totic synchronization phase is given by X^iLi T>iOi{0) / X^ILi + '^sync^ (mod 27r). 
The following remarks concerning Theorem |5.5| are in order. First, notice that 



Theorem 5.5 is in perfect agreement with the results derived in [27 for the case of two 



second-order Kuramoto oscillators. Second, Theorem 5.5 shows that phase synchro 
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nization is independent of the inertial coefficients Mi, tfiereby improving the sufficient 
conditions presented in [T21 Theorems 4.1 and 4.2] and confirming the results in [21[T] 
derived for the infinite-dimensional case. Furthermore, phase synchronization occurs 
almost globally which improves the region of attraction presented in [13 and naturally 
generalizes the result known for the first-order model |43l Theorem 1]. Third, as in 
Section |4j the bound on -ftTcriticai presented in (i) is only tight and may be conservative 
for a particular set of natural frequencies. Since the multi-rate Kuramoto model ( |1.3[ ) 
is an instance of the parameterized system considered in Theorem|5.1[ it has the same 



equilibria and the same stability properties as the scaled Kuramoto model (5.8) (to 
gether with the frequency dynamics (5.9)). Hence, the implicit formulae (3.5)-(|3.6| 



can be applied to the scaled Kuramoto model (5.8) to find the exact critical coupling 
for a given set of natural frequencies. Fourth, we remark that every local bifurcation in 



the multi-rate Kuramoto model ( 1.3 ) is independent of the inertiae Mi since it can be 
analyzed locally by means of the scaled Kuramoto model (5.8). Moreover, the asymp- 



totic magnitude of the order parameter determined by the location of phase-locked 
equilibria is also independent of the inertiae. Fifth and finally, Theorems 5.1 and 5.3 



apply to any variant of the multi-rate Kuramoto model (1.3) that can be written in 



the forced Hamiltonian and gradient form (5.11 with normally hyperbolic equilibria 



For example, the results on almost global phase synchronization for a connected and 
undirected coupling topology [JOl Proposition 3.3.2] and for state-dependent coupling 
weights [41] can be directly applied to the multi-rate Kuramoto model (1.3). 



Based on the results in this section, we conclude that the inertial terms do not 
affect the location and local stability properties of synchronized trajectories in the 
multi-rate Kuramoto model (1.3). However, the inertiae may still affect the transient 



synchronization behavior, for example, the convergence rates, the shape of separatrices 
and basins of attractions, and the qualitative (possibly oscillatory) transient dynamics. 

Proof. We begin by proving the equivalence (i) <4> (ii). By Theorem 5.3 a locally 
exponentially stable synchronized trajectory of the multi-rate Kuramoto model (|1.3|) 



exists if and only if there exists a locally exponentially stable equilibrium of the corre- 



sponding scaled Kuramoto model (5.8). The latter is true if and only if statement (i) 
holds, see Theorem |4.1[ Moreover, Theorem 4.1 asserts that a synchronized solution 
is phase cohesive in A(7t„in). This proves the equivalence (i) (ii). 

We next prove the implication (iv) =^ (iii) . By assumption, there exist constants 



^'sync G IT and 

^sync t 



such that 6i{t) 



"'sync 



' ^sync^ 



(mod 27r), Oiit) 



Wsync, and 



9i{t) = 0, i € {1, . . . , n}. In the phase-synchronized case, the dynamics ( 1.3 ) then read 



as DiUsync = for all I S {1, . . . , n}. Hence, a necessary condition for the existence 



of phase-synchronized solutions is that all ratios uJi/Di 
In order to prove the converse implication (iii) 



cjsync are constant. 

4> (iv), let S = Wsync 

consider the model (1.3) written in a rotating frame with frequency Wsync as 



and 



M'e, + D,0i = sin(6'i -Oj), ie{l,...,m}, 



DiOi = sin(6'j — Oj) , i e {to + 1, . . . , ?i} . 

r) ^- — ^ 



(5.12) 



Note that (5.12) is an unforced and dissipative Hamiltonian system, and the corre 

n ^ 



spending energy function V{9, 6) = ^9'^ AI9 — ^ J2ij=i cos{Oi — 9j) is non-increasing 
along trajectories. Since the sublevel sets of V{9, 9) are compact, the invariance prin- 
ciple |301 Theorem 4.4] implies that every solution converges to set of equilibria. By 
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Theorem 5.3 we conclude that the phase-synchronized equihbrium of (5.121 is locaUy 



exponentially stable if and only if the the phase-synchronized equilibrium of the cor- 



responding scaled Kuramoto model (5.8) with (D^ = is exponentially stable. By [l3l 
Theorem 1], the latter statement is true, all other equilibria are locally unstable, and 
thus the region of attraction is almost global. This concludes the proof of (iii) (iv). 



To obtain the explicit synchronization phase, we sum over all equations ( 5.12[ ) 



to obtain J2iLi Mi9i + X]"=i Di9i = 0. Integration of this equation along phase- 
synchronized solutions yields that X]"=i^i^i(^) = Y^7=i^i^ii^) constant for all 
t > 0, where we already accounted for 9i{t) = for all i G {1, . . . ,n} and all t > 
0. Hence, the synchronization phase is given by a weighted average of the initial 
conditions J^^^i -^i^ii^) / J^i^i -^i- 1^. the original coordinates (non-rotating frame) 
the synchronization phase is then given by -Oi^i(O)/ X]"=i + ^sync^- D 

6. Conclusions. This paper reviewed various bounds on the critical coupling 
strength in the Kuramoto model, formally introduced the powerful concept of phase 
cohesiveness, and presented an explicit and tight bound sufficient for synchronization 
in the Kuramoto model. This bound is necessary and sufhcient for arbitrary distribu- 
tions of the natural frequencies and tight for the particular case, where only implicit 
bounds are known. Furthermore, a general practical stability result as well as various 
performance measures have been derived as a function of the multiplicative gap in 
the bound. Finally, we partially extended these results to the multi-rate Kuramoto 
model and proved that the inertial terms do not affect synchronization conditions. 

In view of the different biological and technological applications of the Kuramoto 
model OEIllinillZllHIZllMllII], similar tight and explicit bounds have to be 
derived for synchronization (as well as splay state stabilization) with arbitrary cou- 
pling topologies, phase and time delays, non-gradient-like dynamics, and possibly 
non-uniform coupling weights depending on state and time. 
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